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ABSTRACT 

We report two epochs of Chandra-ACIS X-ray imaging spectroscopy of the 

nearby bright Type Iln supernova SN 2010jl, taken around 2 months and then a 

year after the explosion. The majority of the X-ray emission in both the spectra 

^ , is characterized by a high temperature (> 10 keV) and is likely to be from the 



l> 



pq 



^ 



^ 



forward shocked region resulting from circumstellar interaction. The absorption 



Q-i. column density in the first spectrum is high (~ 10^^ cm^^), more than 3 orders 

O , of magnitude higher than the Galactic absorption column, and we attribute it to 

^ '. absorption by circumstellar matter. In the second epoch observation, the column 

^ ' density has decreased by a factor of 3, as expected for shock propagation in the 

circumstellar medium. The unabsorbed 0.2 — 10 keV luminosity at both epochs is 

^ ! ~ 7 X 10^^ erg s~^. The 6.4 keV Fe line clearly present in the first spectrum is not 

detected in the second spectrum. The strength of the fluorescent line is roughly 

*sO ! that expected for the column density of circumstellar gas, provided the Fe is not 

highly ionized. There is also evidence for an absorbed power law component 

(^ ! in both the spectra, which we attribute to a background ultraluminous X-ray 

^ ' source. 

L^ . Subject headings: Supernovae: General — Supernovae: Individual (SN 2010jl) — 

;_( ! hydrodynamics — circumstellar matter — X-rays: general 
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Introduction 



Supernova (SN) 2010 il was discovered on 2010 Nov 3.5 (UT) at a magnitude of 13.5 
( JNewton &: PuckettI l2010l ). and brightened to magnitude 12.9 over the next day, showing 
that it was discovered at ai i early phase. P re-discovery observations indicate an explosion 
date in ea rly October 2010 (IStoU et al.ll201ll ). Spectra on 2010 Nov 5 show that it is a Type 
Iln event ( Benetti et al.l 120101 ). The apparent magnitude is the brightest for a Type lln SN 
since SN 1998S. SN 2010jl is associated with the galaxy UGC 5189A at a distance of 50 
Mpc {z = 0.011), implying that SN 2010jl belongs to the class of luminous SNe Iln with 
an absolute magnitude < —20. Hubble Space Telescope (HST) images of the site of the SN 
taken a decade before the SN indicate that, unless there is a chance coi ncidence of a brigh t 
star with the SN site, the progenitor star had an initial mass > 30 Mq (jSmith et al.l 1201 ll ). 
Optical spectra give evidence for a dense circum stellar medium (GSM) expanding arou nd the 
progenitor star with speeds of 40 — 120 km s~^ ( ISmith et al.l 1201 ll ). IStoU et al.l ( 120111 ) found 
that the host galaxy is of low metallicity, supporting the emerging trend that luminous SNe 
occur in low metallicity environments. They determine the metallicity Z of the SN region 
to be < 0.3 Zq. 

Spitzer observations showed a significant infrared (IR) excess in SN 20 10ih indicating 



eithe r new dust format i on or the heating of circumstellar dust in an IR echo (lAndrews et al. 



201ll ). lAndrews et al.l (120111 ) attributed the IR excess to pre-existing dust and inferred a 
massive CS M around SN 2010 il, possibly suggesting a luminous blue variable (LBV) like 



progenitor. ISmith et al.l (120 12[ ) found signatures of new dust formation in the post-shock 



shell of SN 2010J1 from their multi-epoch spectra. While a significant IR excess is present, 
the SN does not show large reddening, i ndicating that t he du st does not have a spherically 
symmetric distribution about the SN (lAndrews et al.l 1201 ll ). The S wift on-board X-ra y 
Telescope (XRT) detected X-rays from SN 2010J1 on 2010 Nov 5.0-5.8 Jimmler et al.lboioh . 
Assu ming a temper a ture o f 10 keV and a Galactic absorption column of Nh = 3.0 x 10^° 



cm 

s"^ in the 0.2 



Immler et al.l (120101 ) obtained an unabsorbed X-ray luminosity of 3.6 ± 0.5 x 10^° erg 



10 keV band. 



After the detection of SN 2010jl with the Swift-XKT , we triggered Chandra observa- 
tions of the SN at two epochs, in Dec 2010 and Oct 2011, and we present the results here 
(Section 12]). We discuss the significant changes in the two Chandra observations taken 10 
months apart in Section 121 
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2. Observations and Analysis 

2.1. Observations 

The Swift detection of SN 2010jl allowed us to trigger our approved Chandra Cycle 11 
program in Dec 2010. We again observed SN 2010jl in Oct 2011 under Cycle 13 of Chandra. 
The first observation (Figure [1]) took place under the proposal # 11500430 starting 2010 
Dec 7 at 04:22:53 Hrs (UT) for an exposure of 19.05 ks and then on 2010 Dec 8 at 00:50:20 
Hrs (UT) for a 21.05 ks exposure. The observations were taken with the ACIS-S without 
grating in a VFAINT mode. A total of 39.58 ks exposure time was used in the data analysis 
and 468 counts were obtained with a count rate of (1.13 ± 0.05) x 10~^ cts s~^. The second 
set of observations (Figure [2]) took place under our proposal # 13500593 starting on 2011 
Oct 17 at 20:25:09 Hrs (UT) for an exposure of 41.04 ks. The observations were again taken 
with the ACIS-S in the VFAINT mode with grating NONE. In a total of 40.51 ks usable 
exposure time, we obtained 1342 total counts, i.e. a count rate of (3.29 ± 0.09) x 10^^ cts 
s~^. We extracted the spectra using CIAO software! and used HEAsofcl to carry out the 
spectral analysis. 



2.2. Spectral Analysis 

As can be seen in Figures [T] and [H the normalized count rate is higher in the Oct 2011 
spectra. This does not necessarily indicate higher intrinsic emission from the SN at the 
later time, because the count rates are absorbed count rates and the unabsorbed emission 
depends on the intervening column density, which may change. In the 5.5 — 7.5 keV range, 
the fluorescent 6.4 keV Fe line is present in the first spectrum but not the second (Figure [3]). 
Here we carry out a detailed spectral analysis of both spectra and determine the significance 
of various emission components along with the Fe 6.4 keV line. 



2.2.1. December 2010 spectrum 

For the December 2010 spectrum, we first fit the high temperature component in the 
2 — 10 keV energy range, where most of the emission in the Chandra energy band lies. We fit 



the spectrum in this range with an absorbed Mekal model (JMewe et al.lll985l : iLiedahl et al. 



http : //asc . harvard . edu/ciao/ 
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19951 ) with metallicity Z = 0.3 Zq (jStoU et al.ll201l[ ). The preferred Mekal temperature 



always hit the upper bound of the temperature allowed in the Mekal model, i.e. 79.9 keV. 
The column density is also very high, with Nh ~ 10^^ cm~^. Since the temperature in the 
best fit models seems high, we checked for the possibility of non-thermal X-ray emission and 
fit a power law model. The column density in this fit is consistent with that of the Mekal 
model; however the photon index is too small to be physically plausible (T = 0.33), and we 
disfavor the non-thermal model. 

When we plot the confidence contours of Nh versus kT, the column density in our fits is 
well constrained, but the upper bound of the temperature is not constrained (Figure H]). We 
established a lower bound on the temperature by assuming a value and finding the goodness 
of fit; T = 8 keV gives a good fit with acceptable x^ value, but not lower values. Thus 8 keV 
is a lower limit on the temperature of the main X-ray emission component of the SN. For 
the lower temperature component in the 0.2 — 2 keV range, we fix the absorption column to 
the Galactic value of 3.0 x 10^° cm~^ since the absorption column for this component is very 
poorly constrained. This component is best fit with a temperature of ~ 2 keV or a power 
law with r = 1.76. 

Figure [1] shows the best fit to the Chandra spectrum taken between 2010 Dec 07.18 
UT and 2010 Dec 08.03 UT, and the best fits are tabulated in Table [H Our complete 
model is the Absorption x Powerlaw + Absorption x(Mekal-fGaussian). 97.9% of the total 
unabsorbed fiux in the Chandra band is carried by the high temperature component. The 
Fe 6.4 keV line carries 1.7% of the total fiux while the \ow-Nh component has only 0.4% of 
the total fiux. The rest energy of the Fe line is 6.39 ± 0.06 keV, which is consistent with the 
Ka line. The equivalent width of the line is EWpe = 0.2 ± 0.1 keV. 



2.2.2. October 2011 Spectrum 

The Chandra spectrum taken between 2011 Oct 17.85 — 18.33 UT is different from the 
first epoch spectrum. We first fit the high temperature component between 2 — 10 keV. The 
temperature in this case also reaches the Mekal model upper limit of 79.9 keV. To consider 
the possibility of non-thermal emission, we fit a power law to the spectrum. However, the 
power law fit to this component yields a photon index of F = 0.45, which is implausible. 

In this spectrum, the column density has decreased by a factor of three; the best fit 
column density for Z = 0.3 Zq is ~ 3 x 10^^ cm~^. The Fe 6.4 keV line is not present, but the 
low-temperature component is there. Because of few data points and a large uncertainty in 
the column density in this component, we again fix the column density to the Galactic value. 
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When we fit the low temperature component with a thermal plasma or a power law model, it 
fits with a temperature of ~ 1 — 2 keV or a power law of ~ 1.7. However, another component 
is still required by the data. When we try to fit this component using the same column density 
as that of the main emission component, we find T = 0.11 keV and a very high unabsorbed 
luminosity ~ 10"^^ erg s~^, which is implausible. Thus we try to fit this component with an 
independently varying column density. The column density associated with this component 
is around 1/4 that of the high Nh component column density. This gives a reasonable and 
physically plausible component and indicates that the fiux is this component is 15 — 20% of 
the total emission. Thus, in this spectrum, we have three components: a hi^-Nn high-T 
component, a hi^-Nu low-T component and a Iovi-Nh component. 

Although the preferred temperature again reaches the upper bound allowed by the 
Mekal model, the error determination shows that there is no upper bound and the lower 
bound to the temperature is 12 keV. Our final model thus is Absorption x Power law + 
Absorption X Mekal+ Absorption X Mekal. Figure |2] shows the best fit to the Chandra October 
2011 spectrum. Table [1] lists the models and best fit parameters. In this case, 81.1% of the 
unabsorbed fiux is carried by the high column density component, 18.2% by the lower-A'^^ 
component and 0.7% is carried by the power law component. 



3. Results and Interpretation 

Here we highlight the main differences between the December 2010 and October 2011 
spectra and discuss the best fit models and their implications. The lower limits on the 
temperature for the two spectra are 8 keV and 12 keV, respectively, showing that a hot 
component is present. The column densities of the main X-ray emission component (high- 
Nh component) are high at both epochs. The column densities at the first and second epochs 
are ~ 10^^ cm~^ and 3 x 10^^ cm~^ (for a metallicity oi Z ^ O.SZ©), respectively. These 
are 3000 times and 1000 times higher than the Galactic column density (3 x 10^° cm^^). 
The high value and variability of N^ point to an origin in the circumstellar medium. The 
excess column density to the X-ray emission is not accompanied by high extinction to the 
supernova, showing that the column is probably due to mass loss near the forward shock 
wave where any dust has been evaporated. This is the first time that external circumstellar 
X-ray absorption has been clearly observed in a supernova. 



Assuming 2010 Oct 10 as the date of explosion ([Andrews et al.ll201ll : lPatat et al.ll201ll ). 
the epochs of the two Chandra observations are 59 and 373 days, respectively. The 0.2 — 10 
keV absorbed fiux at the second epoch (1.1 x 10~^^ erg cm~^ s~^) is higher than that at 
the first epoch (6.5 x 10~^^ erg cm~^ s~^), but this is due to the lower absorption column 



-6- 



density at the second epoch. The actual unabsorbed emission from the SN is constant within 
20 — 30%. At the two epochs, the unabsorbed luminosity in the 0.2 — 10 keV band is ~ 
7 X 10^^ er g s~^, placing S N 2010jl among the most luminous X-ray supernovae yet observed. 
Table 1 of llmmlerl (120071 ) shows that the only other SNe with comparable luminosities are 
Type Iln events or gamma-ray burst associated S Ne at early times. The luminosity of 
3.6±0.5 X 10^*^ erg s~^ found by Swift on 2010 Nov 5 (jimmler et al.ll2010l ) is revised to a value 
close to our Chandra result if N^ ~ 10^^ cm~^ is assumed. In the thermal interpretation, 
the shock velocity can be deduced as Vsh = [lQkT/{3fi)Y^'^ = 7700(A;T/80 keV)^/^ ]^^ g-i^ 
where k is Boltzmann's constant and fi is the mean particle weight. A lower limit of 10 keV 
for the temperature puts a lower limit of the 2700 km s~^ on the shock speed. 

In comparing the observed luminosity to a thermal emission model to find the physical 
parameters, we note that our measu rements give the spec tral luminosity, not the total lu- 
minosity. We use equation (3.11) of iFransson et al.l (119961 ) for the luminosity, adjusting to 
an observed photon energy of ~ 10 keV rather than 100 keV; the G aunt factor is increase d 
to 2 — 3. For the preshock column density, we use equation (4.1) of iFransson et al.l (119961 ). 
These expressions allow for a variation of the preshock density oc r~*, where s is a constant. 
The value s = 2 corresponds to a steady wind and is commonly used, but implies stronger 
evolution than we observe in SN 2010jl. If the circumstellar medium around SN 2010jl is 
due to some presupernova eruptive event, deviation from s = 2 is plausible. Another param- 
eter is m, determined by the expansion of the supernova shock R oc f". For the plausible 
value m = 0.8, we find that s = 1.6 give a reasonable representation of the luminosity and 
Nh evolution. The implied value of the mass loss rate M, normalized to R = 10^^ cm, is 
M_3/f^2 ~ 8v^'^, where M_3 = M/(10~^ M© yr"-*^), Vu]2 is the preshock wind velocity in 
units of 100 km s~^, and v^ is the shock velocity in units of 10'^ km s"^ at the first epoch. 

The high temperature implies that we are observing the forward shock region. The 
physical conditions are such that the forward shock front is close to the cooling regime 
(jChevalier &: Irwinll2012l ). In this case, the luminosity of the forward shock is expected to 
dominate that from the reverse shock and the reverse shock emission may be absorbed by a 
cooled shell, which explains the lack of observational evidence for reverse shock emission. 



In modeling the X-ray absorption in SN 2010jl we have assumed that the absorbing 
gas is not fully ionized. If the circumstella r gas is photoionized by the X-ray emission, the 
absorption is reduced (e.g., IFransson! 1 19821 ). Taking an X-ray luminosity of lO'^^ erg s~^ and 
Ms/vu]2 ~ 8 (at r = 10^^ cm), the ionization parameter is C = L/nr"^ ^ 200; a similar value 
is obtained taking nr ~ Nh ~ 10^^ cm~^ and r = 6 x 10^^ cm for the early e poch. This is in 



a reg ime where the CNO elements may be completely ionized, but Fe is not (JHatchett et al. 
19761 ) . The CNO elements absorb radiation at ~ 1 ke V, so there is the possibility of getting 
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enhanced emission around that energy, as is observed in SN 2010jl. We in vestigated this 



possi bihty by running various cases with the CLOUDY photoionization code ( iFerland et aL 
19981 ). It is possible to obtain cases in which there is a peak at ~ 1 keV, but they had too 



httle absorption in the 1.5 — 3 keV range. We thus favor a background source origin for the 
1 keV emission, especially because the emission remains fairly constant over the 2 epochs. 

The value of N^ in Dec 2010 implies that Tg^ ~ 1, where r^s is the electron scattering 
optical depth thr ough the preshock wind. The Ha line at that time showed roughly symmet- 
ric broad wings (jSmith etal.l 120121 ) that are probably best explained by electron scattering 
in the slow moving wind. IChugail (1200 ll ) estimated that the broad features observed in SN 
1998S require Tes ~ (3 — 4). The required optical depth may be several times that ob- 
served along th e line of sight to the X-ray emission, which could be the result of asymmetry. 
Andrews et al.l (120111 ) found that the column density of dust needed for observed infrared 
emission is larger than that on the line of sight to the supernova, although this is at larger 
radii. 

The Dec 2010 spectrum shows a 6.4 keV feature (Figure [3]), which is identified with the 
narrow Ka iron line. Since the 6.4 keV Fe line arises from neutral or low ionized iron (Fel 
to FeXI), it supports our finding that the radiation field is not able to completely ionize the 
circumstellar gas. A simple estimate o f the expected equival ent width of the Fe line {EWpe) 
can be obtained from equation (5) of iKallman et al.l (120041 ): EWpe = 0.3{Z/Zq)N24: keV, 
where N24 is the circumstellar column density in units of 10^^ cm~^ and the line production 
is due to a central X-ray source in a spherical shell. The expression assumes a flux spectrum 
Fe oc -E~^; Fe oc £'~°'^, which is more appropriate to the hot thermal spectrum here, 
increases EWes by 1.2. The prediction for the metallicity in our case in the early spectrum 
is EWfc = 0.1 keV, and the observed value is 0.2 keV. In view of the uncertainties in 
the model and the observations, we consider the agreement to be adequate. At the second 
epoch, Nh is smaller by a factor of 3, so the strength of the Fe line should be correspondingly 
smaller; this is consistent with the nondetection of the line. The problem with this picture is 
that it assumes the Fe is in the low ionization stage s that produce the Ka line; this requires 
an ionization parameter ( ^ 5 (JKallman et al.ll2004j ). which is below the inferred value. One 
possibihty is that the circumstellar gas is clumped, with a density > 40 times the average; 
another is that the K line emission is from dense gas that is not along the line of sight. 

A thermal fit to the low temperature component implies an absorbing column density 
of (1.37 ±8.44) X 10^*^ cm~^, much less than the column to the hot component and consistent 
with the Galactic column density within the errors. This rules out the possibility that the 
cooler X-rays come from slow cloud shocks in the clumpy GSM or from the reverse shocks. 
The component is also present in the second epoch. It could arise from a presupernova mass 



loss event or from an unrelated source in the direction of the supernova. The components are 
best fit with either a thermal component (T ~ 1 — 2 keV) or a power law with F = 1.6 — 1.7. 
The luminosities of this component in the Dec 2010 and Oct 2011 spectra are 3.5 x 10^^ erg 
s~^ and 5 x 10^^ erg s~^, respectively. The luminosity range and the power law index are 
compatible with a background ultraluminous X -ray source (ULX) , which can typically be 



described by an absorbed power law spectrum (jSwartz et al.l 120041 ). Since the error in the 
flux determination is between 20 — 30%, a factor of 1.4 change in the luminosity at the two 
epochs is consistent with a constant flux. Thus we attribute this component to a background 
source, most likely an ULX, which is associat ed with the blue excess emission region seen 



in the pre-SN HST images (jSmith et al.l 1201 ll ). We examined the HEASARC archives for 



useful limits on such a source, but did not find any. 

The December 2010 spectrum has only one temperature component associated with the 
high column density. However, in the October 2011 spectrum, there are two temperature 
components associated with a high column density, one with temperature > 10 keV and 
another with temperature 1.1 keV. The lower temperature component fits with 1/4 the 
column density of the high temperature component. The fact that the component is absent 
at the first epoch suggests that it is related to the supernova emission. We examined the 
possibility that the emission is the result of reduced absorption due to photoionization of 
the absorbing material, in particular, that lighter atoms are ionized but heavier atoms are 
not. However, we were not able to reproduce the observed emission and the source of this 
emission remains uncertain. 

SN 2010J1 is a special Type Iln supernova because we have been able to catch it in X-rays 
early on with as sensitive an instrument as Chandra and trace the early X-ray evolution. We 
observe dramatic changes over two epochs separated by 10 months. For the first time we 
see clear evidence of external GSM absorption in a supernova. We also find that the GSM 
is not fully photoionized by the SN emission, the SN is very luminous in X-rays, and the 
temperature of the emitting gas is > 10 keV. 

We are grateful to the referee for useful comments. Support for this work was provided 
by NASA through Ghandra Awards GO0-211080X and GO2-13082X issued by the Ghandra 
X-ray Observatory Genter, which is operated by the Smithsonian Astrophysical Observatory 
for and on behalf of NASA under contract NAS8-03060. 
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Table 1. Spectral model fits to the SN 2010J1 spectra 



Spectrum 


Model 


2 / 

X 1^ 


Nh 


Parameter 


Abs. 


fiux 


Unabs. flux 




Mekal 

+ Gaussian 


0.84(23) 


9.70±1;«? X 1023 


kT = 79.9_6g 03 
E = 6.32l«:06 


6.55 
2.43 


X 
X 


10-13 
10-14 


2.44 
4.25 


X 10-12 

X 10-14 


Dec 2010 


+ PowerLaw 
Mekal 


0.87(24) 


3.0 X 1020 (fixed) 
10.59tl-2l X 10^^ 


r = 1.68±0;66 
kT = 8.0 (fixed) 


1.21 
5.57 


X 
X 


10-14 
10-13 


1.29 
3.23 


X 10-14 
X 10-12 




+ Gaussian 






V — R 00+0.06 
E - b.32_o 06 


3.02 


X 


10-14 


5.56 


X 10-14 




+ PowerLaw 




3.0 X 1020 (fixed) 


r = 1.63«;™ 


1.18 


X 


10-14 


1.24 


X 10-14 




Mekal 

+ Mekal 


0.92(75) 


8.381^3 X 1022 


kT = 79.9168.55 
kT = 1.05t°i! 


1.04 
3.76 


X 
X 


10-12 
10-14 


2.13 
4.94 


X 10-12 
X 10-13 


Oct 2011 


+ PowerLaw 
Mekal 

+ Mekal 


0.99(78) 


3.0 X 1020 (fixed) 

3.51+?-52 ^ ^g23 

9.06±|2 X 1022 


r = 1.54t0;^3 
kT = 12 (fixed) 
kT= l.lb+lf^ 


2.37 
9.49 
5.34 


X 
X 
X 


10-14 
10-13 
10-14 


2.42 
2.62 
5.87 


X 10-14 
X 10-12 
X 10-13 




+ PowerLaw 




3.0 X 1020 (fixed) 


r - 1 7fi+0.81 
J- — J-.'0_o.84 


1.97 


X 


10-14 


2.15 


X 10-14 



Note. — Here Nu is in cm— 2, E is in keV, and the fluxes are in erg cm— 2 s-l. The fluxes are given in the 
0.2 — 10.0 keV energy range. The absorbed and unabsorbed fluxes are for that particular component in the model. 
The errors in the fluxes arc 20 — 30%. 
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— Best fit Chandra spectrum of SN 2010jl taken in Dec 2010. The spectrum is best 
a high-T, high-A'^H thermal component, a Gaussian at 6.4 keV and a power law with 
index F = 1.76. 
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Fig. 2. — Best fit Chandra spectrum of SN 2010jl taken in October 2011. Tlie plot is for tfie 
model when the column densities of both temperature components (green dashed line and 
the blue dashed line) are allowed to vary independently in the fit. 
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Fig. 3. — Comparison of the December 2010 and October 2011 spectra around the Fe 6.4 keV 
hne energy. The hne is clearly visible at the early time, but not at the later time. 
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Fig. 4. — The plot shows 68% (black), 90% (red) and 99% (green) confidence contours for 
the Nh versus kT confidence contours plot for the December 2010 spectrum. The column 
density is well constrained while the temperature is unconstrained on the high side. 



